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We study the asymptotic solutions of a version of the Balitsky-Kovchegov evolution with discrete 
steps in rapidity. We derive a closed iterative equation in momentum space. We show that it 
possesses traveling-wave solutions and extract their properties. We find no evidence for chaotic 
behaviour due to discretization. 
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CN ■ I. INTRODUCTION 

g. 

Following the discovery of the rapid rise of the proton structure function at HERA ten years ago p|, a great 
effort has been made to theoretically grasp the high-energy limit of QCD. Within the large-jVc limit and in the 
leading-logarithmic approximation, Balitsky and Kovchcgov have proposed an equation [3] combining the linear 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) growth |3| and the saturation effects due to multiple scattering and gluon 
recombination £|. 

This evolution to high rapidities Y = log(s) is expressed in the colour-dipole framework. Actually, one can view 
_^ , the evolution process as a dipole cascade, associated with a classical branching with a probability distribution given 
£NJ ' by the BFKL kernel. These branching processes appear frequently in statistical mechanics and are well suited to be 
ON ', considered with discrete time steps. In the case of the Balitsky-Kovchegov (BK) equation, rapidity plays the role 
of time and its discretization, acting as a veto in rapidity, can e.g. take into account next-to-leading-order effects 

g: Sil _ 

It has been argued recently Q that discretization of the BK equation may lead to chaotic behaviour of the scattering 
tIh ' amplitude at high energy. This discretization is motivated by the fact that fluctuations around the mean-field (BK) 
Oh, solution can modify the scattering amplitude. This effect can be taken into account by discretizing the number of 
emitted gluons [9j . The evolution can then be considered as a discrete process where the emission of a gluon occurs 
over a rapidity interval Ay. Our aim is to study this problem using a discretization scheme corresponding to an 
_C ■ explicit coordinate-dependent branching process (see further eq. @). 

Recently, a discretized version of the BK equation [Iol | has been proposed, based on its holomorphic separability 
• *h , which is a crucial property of the BFKL kernel in the leading-logarithmic approximation. In addition, the BK equation 
has been shown |TI| to admit traveling-wave asymptotic solutions. It means that, when rapidity increases, the dipole 
density T(r, Y) depends on the dipole size r and the rapidity Y only through the variable rQ s (Y) where the so-called 
saturation scale Q 2 S {Y) cx exp(cw c F) is related to the speed v c of the wave. Both the speed and the tail of the front 
T(rQ s (Y)) can be deduced from the BFKL kernel. 

In this paper, we show analytically and numerically that the discretized version of the BK equation leads to the 
formation of traveling waves and we study their properties. We shall show that the velocity of the wave, i.e. the 
exponent of the saturation scale, increases w.r.t. the continuous equation and the slope of the front decreases. 

Finally, in another study of a discretized version of the BK equation [8| , the authors have suggested the apparition 
of chaos in the evolution. We shall discuss that point using our discretization scheme. 

The paper is organised as follows: in section ^ we present the discrete-time formulation of the BK evolution 
equation, we continue with its solutions and, in particular, the traveling- wave solutions in section IhTI We then finish 
by a conclusion and a discussion of our results. 



t on leave from the fundamental theoretical physics group of the University of Liege. 

^ URA 2306, unite de recherche associee au CNRS. 

*Electronic address: marquet@spht.saclay.cea.fr 

^Electronic address: pesch@spht.saclay.cea.fr 

§ Electronic address: gsoyez@spht.saclay.cea.fr 

"'Electronic address: bialasOth.if.uj.edu.pl 



2 



II. DISCRETIZED EVOLUTION EQUATION 

The impact-parameter-independent Balitsky-Kovchegov (BK) equation for the S'-matrix element corresponding to 
the dipole-target scattering amplitude can be written 



dyS(x 01 ,Y) = J d 2 x 2 ^(x , Xl ;x 2 )[5(x 02 ,r)5(x 21 ,r)-5(x 01 ,r)] 



(1) 



where xy = xj — Xj are the transverse coordinates of the dipoles, Xij = |xy|, r — x$i, and 

TSt \ j2 0! S N C 2 

A(x ,xi;x 2 ) d x 2 = -—^ 2 — 2~ " X2 

is the kernel giving the probability for a dipole of coordinates (xo,xi) to split into two child dipoles of coordinates 
(x ,x 2 ) and (x 2 ,Xi). The BFKL kernel K has the property of holomorphic separability, which means that if we 
introduce 

_ (x 02 )i + ^(Xp 2 ) 2 _ _ (X 02 )i - i(x 02 ) 2 

z n , z , 



'01 



1 01 



the kernel factorises into a function of z times a function of z. 

The discretized version of the evolution kernel is fixed to keep this property 1 : 



Ad(x ,xi;x 2 ) d 2 x 2 



1 a„N r . 



N 2^ \x 2 2 x 2 21 



d 2 x 2 



'01 



with 



N 



d 2 x 2 a s N c ( 



x 01 2tt 2 \x 02 x 21 



a s N c r 2 (a)r(l-2a) 
2tt r 2 (l-a)r(2a) 



(2) 



(3) 



where the regulator a, chosen between and 1/2, allows to separate the real and virtual contributions in Q which 
gives for a rapidity-discretized amplitude 



9y5(x i) 



S' r i+i(xoi) — <SVi(xoi) 
AY 



= N Id x 2 A"d(xo,xi;x 2 )5„(xo2)S' ri (x2i) - NS n (x i). 



One thus fixes the rapidity step AY — 1 /N such that 



Y = nAY = — 

N 



with n 6 IN, 

and the evolution equation for the discretized case takes the iterative form 



5„+i(x i) = J d 2 x 2 K d (x ()7 x 1 ;x 2 ) S n (x 02 )S n (x 21 



This evolution corresponds to a branching process with a probability distribution given by the kernel K d 
Let us derive the evolution equation in momentum space. For this sake, we introduce 

5(,)^ = i _ f (aV kl = J^ 2 » lk ^ 2( °- 1) g(x) = _ Jfee ft ^l°-')T(x) 

W- J <p xe ik.x x 2(a-l) J d 2 xe ik.* x 2(a-l) 

where both T( a )(k) and »§( a )(k) are between and 1. Equation J3J) becomes 

(k) = to (k)l 2 or (k) = 22*0 (k) - [f W (k) 



(4) 



(5) 



(6) 



1 K 

To obtain the regularised kernel one replaces ^ °\ 2 by 



and normalises. This amounts to replace K by NK^ in HI . 
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(a) Rapidity step 



(b) Critical slope 



(c) Critical velocity (crosses: 
numerical results) 



FIG. 1: Theoretical predictions for the traveling-wave properties as a function of the regulator a. 

For the purpose of our studies, a closed equation involving only T^ a ' is required. To achieve this, one first invert 
(0) leading to 



T(x) 



r(o) x 2 
T(l - a) 4tt 



k 2 x 2 



f (a) (k) 



T(2a) x 2 
T(l -2a) 4^ 



d 2 ke- lk '- 



k 2 x 2 



T (2a) (k), (7) 



where the second equality is simply obtained by replacing a by 2a. Inserting this into equation Q allows to express 
(k) in terms of T^?\ (k) . The evolution equation JJJJ finally takes the form 



3& ) 1 (k) = 



r(i-o)r(o + |) 2 4q /• d 2 p 



r(o)r(|-o) 2tt y (k- P 



fc 2 (k- P ) 2 



2fW(p)- [fW(p) 



(8) 



It is clear from ||SJ that if T„ (p) depends only on the modulus p of p, T^9 1 (p) also features that property. Thus, 

starting with an angular- independent initial condition Tg a ^(p), this property is conserved throughout the evolution. 
Hence, one can integrate out the angular <p dependence in (JSJ). Introducing the integration variable sin 2 (y?/2), one 
obtains 



lQ r(l-a)r(a+|) f°° pap 



r(a)T(i-a) 7 + 



r fe 2 ( fc+p) 2 

P 4 



2 F, ( 'I,l- a ;l ;77 i^) J 2 fW( P )- [fW^ 2 



(fc + pf 



( 9 ) 

Before considering the solutions of this process, let us show how one can reconstruct the physical amplitude from 
T^ a '. The physical T-matrix in momentum space is defined as usual by 



(10) 



If we combine this with equation Q , one gets 



T n (k) 



pdp 
o (fc+p) 5 



(fc + p) 2 

2 



^(ii-a;i;A)«. 



(fc + p) 2 



(ii) 



Equation together with the transformation (TTJ, constitute a closed iterative description of the discretized 
evolution. 



III. ASYMPTOTIC PROPERTIES: TRAVELING- WAVE SOLUTIONS 



In this section, we discuss the solutions of the branching process © and compare them with the properties of the 
Balitsky-Kovchegov equation . 
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To begin with, let us recall the important properties of the asymptotic solutions of the BK equation. It has been 
observed recently [Tl| that the BK equation admits traveling waves as asymptotic solutions: at large rapidities, the 
amplitude takes the form of a wave of speed v c 

f(k, Y)=T [\og{k 2 ) - av c Y] . 

In physical terms, this means that T depends only on the ratio k/Q s (Y), i.e. satisfies the geometric-scaling property, 
where the saturation scale is given by Q 2 oc exp(av c Y) . 

One can even go beyond this result, valid for infinitely large Y , and compute pre-asymptotic corrections to the 
speed of the wave. Expressed in terms of the saturation scale, one obtains [Tj 



d Y log [Ql(Y)} =av c --Ll + A 1 



2tt 



2 lc Y 2 7c 2 y a X "(lc)Y 3 

where xil) is the eigenvalue of the BFKL kernel and the critical parameters j c and v c are determined by 

I, \ Xilc) 

Vc = X (7c) = • 



In addition, at large rapidities, the tail of the front is given by 



T{k) rZj^y mY) 



k 2 



Q 2 S (Y) 



exp 



1 



2a X "(lc)Y 



log 2 



l(Y) 



(12) 



(13) 



(14) 



Now that we have summarised the properties of the asymptotic solutions of the BK equation, let us consider the 
case of the discretized equation (0J . More precisely, we shall prove that one can also infer traveling- wave solutions in 
this equation. 

In order to achieve this task, we have to show that the solution of the linear equation (T = 1 — S) 

T n+ i(x i) = J d 2 x 2 K d (xo,x 1 ;x 2 ) [T„(x 02 ) + T„(x 2 i)] 
can be expressed as a superposition of waves. Going to the Mellin representation 

J 2l7T 

one easily gets 

r„(7) = r ( 7 )exp[iv x ( a )(7)y 

with 



X (Q) (7) = log 



T(a + 7) T(l - a) r(l - 2a - 7) T(2a) 



r(a)r(i 



7)r(l-2a)r(2a + 7 ) 



(15) 



The kernel X which is the logarithm of the Mellin transform of the kernel Kd (see (j2Jl), describes the linear 

growth, as does the BFKL kernel x(l) i n the continuous case. One can thus predict that equation admits 
traveling waves as asymptotic solution with the critical exponent 7 C and the critical speed v c fixed by equation 1131) 
with Nx^ a Hl) instead of x(t)- The rapidity dependence of the saturation scale and the tail of the front are given by 
equations H12J) and (|14(l . where the critical parameters are now obtained from the kernel Nx^(l) : 



(16) 



In figure ^ we represent the main properties of the traveling waves as a function of the regulator a. The value of 
the rapidity step (or, more precisely, aAY) shows a maximum of 0.258 for a rs 1/3. When a goes from to 0.5, the 
asymptotical slope of the wave, given by 7c , decreases 2 while its speed increases. In the limit a — > 0, one recovers the 
parameters from the BK equation: 7 C pa 0.6275 and v c « 4.88. 



2 Up to a small correction of order 10 4 , we observe that 7c(<i) — 7c(0)(l — 2a). 
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log(fc) 



FIG. 2: Traveling-wave pattern from discretized evolution. We have fixed a = 0.3 and a = 0.2, corresponding to rapidity steps 
of AF = 1.264. The amplitude f n (k) is shown, from left to right, for n = 0, 1, . . . , 14. 




FIG. 3: Comparison of the traveling-wave front for different values of the regulator a. We observe that the slope decreases as 
a increases. 
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FIG. 4: Speed of the traveling wave. The solid curve shows the theoretical prediction 1121 for 9ylog(Q^). The gray area 
corresponds to numerical computation using different extractions of the saturation scale (see text, eq. 11710 : the lower(higher) 
border corresponds to T(Q S ) — 0.1 (0.001). The theoretical result for the continuous case has been added for comparison. 
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In order to check those predictions, we have performed a numerical analysis of equation © . We start with a initial 
condition for the physical amplitude given in coordinate space by 

S(x) = exp (-sb") , 

where the parameter p is chosen such that both S(x), T^- a \k) and T(k) are positive. Changing p allows us to discuss 
the effect of the initial condition on the formation and properties of the traveling wave. From this initial condition, 
one can compute T^ a \k) using (jSJ). The evolution is then performed by iterating the branching process © and the 
physical amplitude is recovered through JTTJ. The value of the strong coupling constant has been kept fixed a = 0.2. 

First of all, let us discuss the formation of the traveling wave. To emphasise clearly the effect of discretization, 
we have chosen a = 0.3, a = 0.2 and 3 p = 1.2, leading to a large rapidity step AY = 1.264 (for a = 0.1, this gives 
AF w 0.5, see figure QJ. The result of the evolution is plotted in figure [21 where each curve corresponds to one step 
in rapidity i.e. to To(fc), . . . ,Tn{k). We obviously see from the figure that the trace of the initial condition is lost 
during the evolution and that a traveling wave is formed. One can also show that if the initial condition is not steep 
enough, the slope of the initial condition is preserved by the evolution as for the continuous case. 

If we consider different values of the regulator a, the slope of the resulting traveling wave, given by 7 C provided the 
initial condition is steep enough, shall depend on a. This slope, extracted from JTBJ, is expected to decrease when a 
increases. In figure |3 we show the form of the front obtained after evolution for different values of a, and it clearly 
appears that the slope decreases significantly for large values of a. 

Finally, we have to check that the rapidity dependence of the saturation scale is in agreement with equation l|12|l . 
To extract the saturation scale from the numerical simulations we solve 

f(Q s (Y),Y)=f . (17) 

The level To at which we extract the saturation scale is not fixed, but it only introduces subleading corrections 
to the saturation scale when rapidity increases. In our analysis, we have chosen To between 0.001 and 0.1, which 
is characteristic for the region where the traveling wave is formed. We can then plot the numerical value for the 



saturation scale together with the predictions from 1)12(1 . The result is shown in figures 1(c) and 0] for different values 
of the regulator a and for p = 1.2, together with the theoretical prediction for the continuous BK equation. In addition 
to the perfect agreement between the theoretical expectation and the numerical results, we observe that the speed of 
the traveling wave becomes larger and larger as a increases. 



IV. DISCUSSION AND CONCLUSION 



Throughout this paper, we have seen that the property of holomorphic separability gives a natural way to discretize 
the BK equation in rapidity through the introduction of a regulator a which allows to adjust the rapidity step. The 
evolution then turns into a discrete branching process which can be cast in a closed form © from which we can easily 
recover the physical dipole amplitude l|ll|) . 

We have studied the asymptotic properties of this discretized equation both analytically and numerically. An 
analysis of the new evolution in Mellin space predicts the formation of traveling-wave solutions as it is the case for 
the continuous BK equation. We have tested numerically that these traveling waves indeed appear with the correct 
rapidity dependence of the saturation scale Q 2 S and the correct behaviour of the front. The most important effects 
of discretization is to increase the speed of the wave and to tame the fc 2 -dependence of the tail of the front when a 
increases (a — > corresponds to the continuous limit). 

It is important to mention that, in comparison with the approach proposed recently in Ref. |8(, our evolution 
equation does not generate chaos. We do believe that this is due to a more careful treatment of the rapidity steps. 
Indeed, if one discretizes the Y-derivative of the BK equation using steps in ctY of size 5 and keep only the lowest 
order expression for the BFKL kernel x(t) ~ x(l/2), one gets 

T n+1 = [l + dx(l/2)]T n -Tl 

Defining U = [1 + 5x(l/2)] -1 T, the evolution takes the form of the famous logistic map ^3 

U n+X = XU n (l - U n ), 



3 The choice of the initial condition is due to physical constraints, imposing a finite range. For example, for p = 2 the Fourier-transformed 
amplitude T^ a '(k) is not positive definite, while p = 1 corresponds to a subcritical slope which is not the case in QCD. 
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with the Malthusian parameter A = 1 + 5^(1/2). In 8], the authors adopt 5 = 1, corresponding to large rapidity 
steps AY = 1/a « 5, giving A = 3.77. This value being larger than the accumulation point A ap « 3.56994567, the 
evolution generates chaos. In the approach developed in this paper, the value of 5 is given by l/N which depends on 
a. One has noticed that the maximal value of the step size is 5 ks 0.258, reached for a w 0.333. This corresponds to 
a Malthusian parameter of 1.715. Since this is smaller than 3, we lay in the 1-cycle region and no chaos is observed. 
Note finally that our evolution equation , which does not require any approximate expansion of the kernel, is itself 
very similar to the logistic map up to a specific convolution (see formula ©)• 
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